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ABSTRACT 



Aims. To compare two different diagnostics for estimating stellar masses in early-type galaxies and to establish their level of reliability. 
In particular, we consider the well-studied sample of 15 field elliptical galaxies selected from the Sloan Lens ACS (SLACS) Survey 
(z = 0.06 - 0.33). We examine here the correlation between the stellar mass values, enclosed inside the Einstein radius (/?Ein) of each 
lens, based on analyses of lensing and stellar dynamics combined and based on multiwavelength photometry spectral template fitting. 
Methods. The lensing+dynamics stellar mass A/]* n+d (< Re™) is obtained from the published SLACS Survey results, assuming a 
two-component density distribution model and a prior from the fundamental plane on the mass-to-light ratio for the lens galaxies. The 
photometric stellar mass Af* hot (< /?Ein) is measured by fitting the observed spectral energy distribution of the galaxies (from the SDSS 
multi-band photometry over 354 - 913 nm) with composite stellar population templates, under the assumption that light traces stellar 
mass. 

Results. The two methods are completely independent. They rely on several different assumptions, and so, in principle, both can 
have significant biases. Based on our sample of massive galaxies (^°S^ hot (< /?h„) - [10.3, 11.5]), we find consistency between the 
values of MJ. n+dyn (< /?Ein) and M^ bot (< /?Ein)- We obtain a Pearson linear correlation coefficient of 0.94 and a median value of the 
ratio between the former and the latter mass measurements of 1.1 ± 0.1. This suggests that both methods can separately yield reliable 
stellar masses of early-type galaxies, and confirms that photometric mass estimates are accurate, as long as optical/near-IR rest frame 
photometry is available. 
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1. Introduction 

An estimate of the stellar mass component in galaxies is inter- 
esting for several different reasons. In detail, by combining or 
comparing photometric stellar mass estimates, obtained by spec- 
tral energy distribution (SED) fitting methods, with dynamical 
or lensing measurements, it is possible to study the radial dis- 
tribution of dark matter (e.g., Ferreras et al. I2005I Napolitano 
et al. |2005l l, to investigate the relationship between stellar and 
total mass (e.g., Lintott et al. 2006; Rettura et al. I2006l l. and to 
test hierarchical structure formation models (e.g., Nagamine et 
al. |2004| De Lucia et al. 2006 ). Interestingly, Treu & Koopmans 
(2004 ) have proved that the stellar mass fraction in elliptical lens 
galaxies can also be estimated with a joint lensing and dynamical 
analysis. 

Although it is common to measure stellar masses through 
these techniques, only a few studies have been performed to 
check the reliability of each method (e.g., Drory et al. [2004). 
Further investigations are therefore important to probe the con- 
sistency of the different techniques. 

Throughout this work we assume Hq - 70 km s _1 Mpc~', 
Q m = 0.3, and Q A = 0.7. 



2. The SLACS sample 

In this Letter, we focus on a uniformly selected sample of 15 
massive field early-type galaxies taken from the SLACS Survey 
(for more details on the selection procedure, see Bolton et al. 
2006). Table Q] summarizes the relevant photometric and spec- 
troscopic properties of the galaxy sample. The lens galaxies have 
a redshift between 0.06 and 0.33, HST F435W and F814W im- 
ages, usriz magnitudes and stellar velocity dispersions from the 
SDSSQ- They are luminous red galaxies (LRG; Eisenstein et al. 
120011 ) with properties similar to those of non-lensing early-type 
galaxies: redshifts, stellar velocity dispersions, stellar popula- 
tions, and mass density profiles (see Treu et al. 12006; Koopmans 
et al. 120061 



3. Measuring stellar masses 

Here, we describe how the stellar mass of the galaxies of our 
sample is measured using two independent diagnostics. 
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3.1. Lensing+dynamical estimates 

Strong gravitational lensing provides the most accurate estimate 
of the total (stellar +dark) projected mass of a lens galaxy in- 
side the Einstein radius. It has been shown by Treu & Koopmans 
(2004) that by combining lensing measurements with spatially 
resolved kinematic profiles in elliptical galaxies, the stellar and 
dark components can be separated precisely. If the velocity dis- 
persion of stars is known only from a single (fiber) aperture, 
some information on the stellar mass fraction (/„) inside /?Ein can 
still be obtained. This particular analysis has been performed on 
the SLACS sample by Koopmans et al. (12006b . The results are 
shown in Table [2] We summarize here the main steps and as- 
sumptions: 

- The mass distribution of each lens galaxy is described in 
terms of a two-component spherical and isotropic model. 
It consists of a Hernquist density profile, scaled by a stel- 
lar mass-to-light ratio M*/L, for the stellar component, and 
a power-law density profile, for the dark matter component 
(Pd = Pd,o r~ y ). 

- The lensing measurement of the total projected mass en- 
closed inside the Einstein circle, M^ +d ^ n (< Reir), is used to 
eliminate pa,o- Thus, for any given set {M*/L, y), the spher- 
ical Jeans equation can be solved to determine the line-of- 
sight stellar velocity dispersion as a function of radius. 

- A likelihood function is defined by comparing the model pre- 
dicted with the observed aperture stellar velocity dispersion. 
A prior on the stellar mass-to-light ratio, based on the funda- 
mental plane (in the B band) and corrected for passive evo- 
lution, is considered before marginalizing on the two free 
parameters of the model (M*/L and y). 

- The stellar mass fraction is calculated as the ratio between 
the maximum likelihood and the maximum allowed values of 
M*/L. This latter quantity is obtained under the assumption 



that ^en+dyn^ *Ein) 



KU„(^Ei„). 



In principle, this method can give values of /» bigger than 1, as 
can be seen in a few cases in Table|2] Nevertheless, these values 
are consistent with 1, so that the previous assumptions are not 
in doubt. Finally, Koopmans et al. (2006) point out that more 
realistic values of M*/L should take into account the dependence 
on the mass of the galaxy. 

The stellar mass obtained from the combination of gravi- 
tational lensing and stellar dynamics measurements inside the 
Einstein radius, M* en+d (< ^EinX is then estimated by multiply- 
ing the total mass by the stellar mass fraction. This is shown in 
the following equation 



M, 



len+dyn*^ ^Ein) = /. X M^ +dyn (< ^Ein) , 



(1) 



and the results for the SLACS sample (Koopmans et al. 2006) 
are reported in Table [2] 

3.2. Photometric estimates 

The flux of our galaxies has been measured in five different 
bands (see Table Q]), from 354 to 913 nm (ugriz filters of the 
SDSS). In order to obtain unbiased galaxy colors, in the absence 
of color gradients, a de Vaucouleurs profile, 



I(R) = 7 exp 



- 7 < 
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Fig. 1. SED and best-fit model of the lens galaxy SDSS 
J0912+0029 at z = 0.1642. The circles with the error bars repre- 
sent, from left to right, the observed total flux densities measured 
in the u, g, r, i, and z SDSS passbands. On the top, the best-fit 
values of the age (T), the characteristic time of the SFH (t), and 
the mass (M* hot ) of the galaxy are given. On the bottom right, 
the inset shows the 68% and 99% confidence regions for T and 

^^phot ' 



(R e being the standard optical effective radius), is fitted in the r 
band of each object and this model is then varied, only in ampli- 
tude, for the other bands (after convolution with the correspond- 
ing point spread function). The resulting magnitudes are called 
modelMag in the publicly available SDSS catalog and corre- 
spond to magnitudes measured through equivalent apertures in 
all bands. 

We derive the total photometric stellar mass (M* not ) by fit- 
ting the observed SED with a three-parameter grid of composite 
stellar population (CSP) models. We use Bruzual & Chariot's 
(2003 ) templates at solar metallicity, assuming a Salpeter ( 1955 ) 
time-independent initial mass function (IMF) and a delayed ex- 
ponential star formation history (SFH) parametrized by a charac- 
teristic timescale (t). The other two free parameters are the age 
of the model (T), which is constrained by the age of the Universe 
at the galaxy redshift, and the stellar mass. The uncertainties on 
the best-fit parameters are estimated by projecting the joint prob- 
ability density distribution onto the corresponding axes. As an 
example, Fig. [T] shows the SED of a galaxy in our sample and 
the best-fit model. The estimated stellar mass values are given 
in Table [2] We emphasize that the accuracy and homogeneity of 
the SDSS photometry, particularly the lack of significant system- 
atics due to calibration or aperture corrections, ensure accurate 
photometric mass estimates. 

The photometric stellar mass inside the Einstein radius, 
M* hot (< /?Ein), of Table|2]is inferred by multiplying the total mass 
by an aperture factor (/ ap ): 



M; hot (< R Em ) = /a, 



x M. 



phot ' 



(3) 



This last factor represents the fraction of light enclosed inside 
the Einstein radius, with respect to the total light of the galaxy 
parametrized by the de Vaucouleurs profile, and is explicitly 
given by the following expression: 



(2) /a 



ap 



J^ Ein I(R)RdR 
f™I(R)RdR 



(4) 
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The assumption implicit in using Eq. is that the stellar mass 
is traced by the light distribution. 

Finally, we note that stellar mass estimates depend on the 
adopted IMF and SFH, but not on the different stellar population 
models (e.g. Bruzual & Chariot [20031 vs. Maraston 12505] ) . We 
will come back to these points in the next section. 

4. Comments and conclusions 

In Fig. [2] we compare our photometric stellar mass estimates 
with the lensing+dynamics ones of Koopmans et al. (2006 ). This 
plot shows that the two different mass measurements are statis- 
tically correlated and consistent within the error bars. In par- 
ticular, the value of the Pearson linear correlation coefficient is 
p = 0.94, and the best-fit correlation line yields: 



io 2 - 30±168 xM; hot (< * Ein )™- 15 . 



(5) 



The median value of the ratio (q) between the lens- 
ing+dynamical and photometric stellar mass estimates is con- 
sistent with unity (1.1 + 0.1). The ratio q does not show any cor- 
relation with galaxy colors, hence excluding a possible source of 
systematic errors in the photometric mass estimates. 

No significant difference in the relation between M* en+dyn (< 
/?Ein) and M* hot (< ^Ein) is observed deriving photometric stel- 



lar masses with Maraston's (2005 ) CSP models, as the two stel- 
lar population models differ remarkably in a near-IR regime at 
wavelengths longer than the ones probed by the SDSS filters. 
This result is in agreement with that reported by Rettura et al. 
(120061 ). By choosing Kroupa J200TT ) or Chabrier's J2003]) IMFs, 
the photometric mass measurements are lowered in such a way 
that the slope of the best fit is consistent with that of Fig. |2l 
but the average ratio between the two mass estimates is consid- 
erably smaller than one. Moreover, we note that the effect of 
contamination caused by the lensed objects in the measurements 
of the fluxes of a lens galaxy is small. This is usually more rel- 
evant in the bluer filters, which are known to be less sensitive 
to photometric mass estimates, and the measurements of total 
magnitudes through de Vaucouleurs profile fitting should further 
reduce this source of uncertainty. 

Values of q slightly larger than one, like those observed and 
shown in the inset of Fig. |2j may be explained by possible un- 
derestimates of A^* hot (< ^EinX which can be ascribed to two dif- 
ferent phenomena occurring in the galaxies: dust extinction and 
metallicity values lower than the solar one. In detail, both effects 
tend to give lower IR fluxes, which then result in lower mass es- 
timates. Nevertheless, several tests have supported the validity 
of the dust-free and solar metallicity model assumptions (e.g., 
see Rettura et al. l2006b . 

Finally, despite a number of assumptions, we conclude that 
the good agreement between the two mass estimators, within 
their respective uncertainties, is a very reassuring result. This 
makes the presence of strong biases in one of the two meth- 
ods very unlikely, and allows us to use either of them indepen- 
dently to measure reliably stellar masses. Although this study is 
based on a low redshift sample, we expect photometric mass es- 
timates to be also accurate to high redshift, as long as the same 
optical/near-IR rest frame bands are covered. 
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Fig. 2. Comparison of the lensing+dynamical and photomet- 
ric stellar masses (measured inside the Einstein radii) for the 
SLACS sample of gravitational lens early-type galaxies. The 
plot shows M* n+dyn (< R E in) vs. M* hot (< R E i n ), with the best 
fit correlation (solid) and the M;* en+dyn (< R Eia ) = M* hot (< ^Ein) 
(dotted) lines. The diamonds indicate upper limits for the lens- 
ing+dynamical stellar masses obtained by assuming that the 
measured masses enclosed inside the Einstein circles are only 
stellar. In the inset, we plot the ratio (q) between the lens- 
ing+dynamical and photometric stellar mass estimates as a func- 
tion of M* h J< R Ein ). 
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Table 1. The relevant spectroscopic and photometric measurements of the SLACS sample. 
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Notes - Magnitudes are extinction-corrected modelMag (AB) from the SDSS. 
References - Treu et al. ( 20061; Koopmans et al. d2006b . 



Table 2. The lensing+dynamical and photometric mass measurements of the SLACS sample. 
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